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INERTIAL  RANGES  IN  TWO-DIMENSIONAL  TURBULENCE 

Robert  H.  Kraichnan 
Peterborough,  New  Hampshire 


Two-dimensional  turbulence  has  both  kinetic  energy  and  mean-square 
vorticity  as  inviscid  constants  of  notion.  Consequently  it  admits  two 
formal  inertial  ranges,  E(k)  ^  and  F(k)  n.  3,  where  r  is 

the  rate  of  cascade  of  kinetic  energy  per  unit  mass,  n  is  the  rate  of 

cascade  of  mean-square  vorticity,  and  the  kinetic  energy  per  unit  mass 

« 

f  m 

is  Jq  E(k)dk.  The  -S/3  range  is  found  to  entail  backward  energy  cascade, 

from  higher  to  lower  wavenumbers  k,  together  with  zero  vorticity  flow. 

* :.e  -3  range  gives  an  upward  vorticity  flow  and  zero  energy  flow.  Tha 

paradox  in  these  results  is  resoVed  by  the  irreducibly  triangular  nature 

of  th9  elementary  wavenumber-interactions.  The  formal  -3  range  gives  a 

nonlocal  cascade  and  consequently  must  be  modified  by  logarithmic  factors. 

If  energy  is  fed  in  at  a  constant  rate  to  a  band  of  wavenumbers  k,  and 

the  Reynolds  number  is  large,  it  is  conjectured  that  a  quasisteady  state 

results  with  a  -5/3  ra'-.ge  for  k  <<  k.  and  a  -3  range  for  k  *>  k.,  up  to 

X  with  time 

the  viscous  cut-off.  The  total  kinetic  energy  increases  steadilyAas  the 
-5/3  range  pushes  to  ever-lower  k,  until  scales  the  sire  cf  the  entire 
fluid  are  strongly  excited.  The  rate  of  energy  dissipation  by  viscosity 
decreases  to  zero  if  kinematic  viscosity  is  decreased  to  zero  with  other 

parameters  unchanged. 
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1.  INTRODUCTION 

The  vorticity  of  each  fluid  element  is  an  inviscid  constant  of 
notion  in  two-dimensional  incompressible  flow.  Therefore  the  mean-square 
vorticity  as  well  as  the  kinetic  energy  per  unit  mass  art  inviscid  con¬ 
stants  in  two-dimensional  isotropic  turbulent.  A  number  of  authors 
have  studied  twe-din«nsional  turbulence  theoretically,  and  it  is  recognized 

that  the  vorticity  constraint  has  profound  effects  on  inertial  energy 
1-9 

transfer.  In  contrast  to  the  predominantly  one-way  flow  of  energy 
familiar  in  three  dimensions,  .ransfer  upward  in  wavenumber  must  be 
accompanied  by  comparable  or  greater  downwat-d  transfer. 

According  to  the  two-dimensional  Navier-Stokes  equation,  the  inter¬ 
action  of  each  triad  of  wavenumbers  k,  p,  q  individually  conserves  both 
energy  and  squared  vorticity.  In  order  to  separate  off  questions  involving 
the  localnes3  of  energy  transfer,  suppose  that  all  triad  interactions  for 
which  the  smallest  of  the  three  wavenumbers  is  less  than,  say,  one-half 
the  largest  are  aroitrarily  eliminated  from  the  dynamical  equations. 

Divide  the  wavenumber  range  into  half-octave  segments,  so  that  all  the 
triad  interactions  left  in  the  equations  either  connect  nearest -neighbor 
segments  or  are  contained  entirely  within  a  single  segment.  Fig.  1-a 
illustrates  the  segmentation  and  the  way  in  which  the  triad  interactions 
connect  nearest  neighbors.  In  order  for  both  energy  and  squared  vorticity 
to  be  conserved,  the  net  transfer  by  each  triad  interaction  must  either  be 
out  of  the  middle  wavenumbei  into  both  smallest  and  largest  wavenumbers, 
or  vice  versa.  The  arrows  in  Fig.  1-a  are  arbitrarily  drawn  for  the  case 
of  outflow  from  the  middle  wavenumber. 
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In  three  dimensions,  the  triad  interactions  can  be  likened  with  some 
success  to  pair  interactions  (Fig.  1-b).  This  is  inadmissible  in  two  di¬ 
mensions  because  pair  interactions  cannot  transfer  both  energy  and  squared 
vorticity  conservatively  between  unequal  wavenumbers. 

There  is  no  simple  general  relation  between  energy  and  squared- 
vorticity  transfer.  Segment  n  in  Fig.  1-a  is  connected  to  the  lower-wave- 
nunber  segment  n-1  by  two  kinds  of  triad  interactions:  those  with  a  pair 
of  wavenumbers  in  n  and  those  with  a  pair  of  wavenumbers  in  n-1.  The 

former  transfer  squared  vorticity  and  kinetic  energy  per  unit  mass  from 

2 

n-1  to  n  in  ratios  <  2k  ,  where  k  is  the  wavenumber  which  separates  the 

2 

segments.  The  latter  transfer  these  quantities  in  ratios  >  2k^  .  The 
net  rates  of  energy  transfer  per  unit  mass  c  and  squarcd-vorticity  transfer 
n  from  below  k^  to  above  k^  depend  on  the  sig..o  and  relative  strengths  of 
the  two  kinds  of  interactions.  For  example,  if  the  interactions  having 
a  single  wavenumber  in  n  are  sufficiently  strong  compared  to  those  having 
?  single  wavenumber  in  n-1,  it  is  possible  for n to  be  positive  while 
simultaneously  e  is  negative. 

An  important  inference  can  be  made  for  similarity  cascades  where  a 

k-independent  total  contribution  to  c  is  made  by  all  triads  whose  ratio  of 

largest  to  smallest  wavenumber  falls  below  some  arbitrary  limit.  The  triple 

moments  can  be  chosen  to  construct  such  ranges  at  a  given  instant.  Whether 

they  are  self -preserving  is  another  matter.  By  similarity,  the  vorticity 

2 

cascade  rate  must  have  the  form  n  =  2Ak  c  with  A  k-independent.  But  the 
rate  itself  must  be  independent  of  k.  Otherwise,  the  outflow  of  squared 
vorticity  from  each  segment  would  not  equal  the  inflow,  which  would  violate 
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vorticity  conservation  since  with  k-indepcndent  c  the  rate-of-change  of 
excitation  intensity  is  instantaneously  zero  at  each  k.  The  only  possible 
resolution  is  A  :  0,  That  ist  the  rate  of  squared- vorticity  cascade  is 
identically  zero  in  a  similarity  cascade  where  c  is  independent  of  k. 

This  is  corroborated  by  formal  analysis  in  Sec.  2. 

The  roles  of  energy  and  squared  vorticity  are  interchangeable  in 
the  preceding  argument.  If  there  is  a  similarity  range  with  k-independent 
contribution  to  n  from  the  triad  interactions  whose  ratio  of  largest  to 
smallest  wavenumber  rails  within  some  limit,  then  the  contribution  of 
those  triads  to  c  ..us*  be  identically  zero  within  the  range.  This  means 
that  two  kinds  of  putative  inertial-transfer  similarity  ranges  must  be 
investigated:  energy-transfer  and  vort icity-transfer  ranges.  Kolmogorov’s 
assumption  that  the  eneigy  spectrum  E(k)  depends  only  on  k  and  c  1-ads  to 

E(k)  *  Cc2/3k‘5/3  (1.1) 

in  two  dimensions  as  well  as  in  three.  Here  E(k)  is  defined  so  that  t  e 

mean  kinetic  energy  per  unit  mass  is  /Q  E(k)dkvand  C  is  a  constant  whose 

value  can  depend  on  the  dimensionality.  The  alternate  assumption  that  the 

on 

squared-vorticity  spectrum  2k‘E(k)  dependsAonly  n  and  k  yields 

E(k)  *  C»n2/3k‘3,  (1.2) 

where  C'  is  another  constant.  In  Sec.  2  it  is  shown  that  (1.1)  and  (1.2) 
each  satisfy  both  conservation  laws.  Necessary  conditions  for  the  physical 
realizability  of  these  similarity  ranges  are  that  the  transfer  processes 
be  sufficiently  local  in  wavenumber  when  all  triads  are  admitted.  This  is 
discussed  in  Sec.  where  it  is  noted  that  the  -3  range  fails  by  logarithmic 
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factors  to  be  suff'ciently  local  and  therefore  must  be  modified  by  factors 
with  logarithmic  k  dependence. 

Both  experiment  and  general  statistical-mechanical  considerations 
indicate  that  the  energy  cascade  through  the  inertial  range  is  frtro  lower 
to  higher  wavenumbers  in  three-dimensional  turbulence.  The  intensity  at 
high  wavenumbers  is  suppressed  by  viscosity  and  it  is  natural  to  expect 
a  net  transfer  toward  these  wavenumbers  from  the  strongly  excited  low 
wavenumbers,  .n  two  dimensions,  the  vorticity  constraint  drastically 
changes  matters.  A  given  triad  interaction  spreads  the  excitation  in 
wavenumber  space  if  it  gives  a  net  flow  out  of  the  middle  wavenumber  into 
the  small  and  large  wavenumbers.  The  reverse  flow  concentrates  the  exci¬ 
tation.  Spreading  of  the  excitotion  by  the  triad  interactions  would 
seem  to  be  the  more  plausible  state  of  affairs.  Some  supporting  evidence 
is  presented  in  Sec.  3.  If  the  triad  interactions  do  spread  the  excitation 
in  wavenumber  space,  then  it  is  proved  in  Sec.  3  that  the  -S/3  range  yields 
c  <  Oj  that  is,  the  ener^  cascade  is  downward  in  wavenumber.  The  -3 
range  under  the  same  condition  yields  vorticity  cascade  upward  in  wave- 
number  (n  >  0). 

If  the  directions  of  cascade  are  as  just  described,  the  -5/3  range 
could  serve  to  remove  energy  from  an  input  range  of  wavenumbers  down  toward 
zero  wavenumber,  while  the  -3  range  could  carry  vorticity  up  to  the  dissi¬ 
pation  range.  Thus  both  ranges  could  exist  simultaneously.  This  conjec¬ 
ture  is  discussed  in  Sec.  4. 
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2.  FORMAL  ANALYSIS  OF  THE  .SIMILARITY  RANGES 

Let  the  flow  be  confined  in  a  cyclic  box  of  side  D  and  expand  the 
velocity  field  in  Fourier  series  so  that  the  incompressible  Navier-Stokes 
equation  becomes 

0/3t  ♦  vk2)u,(k)  =  -ik  (6.  ,-k,k./k?)  1  .  u.(p)u  (q),  (2.1) 

1  ^  Hi  1]  1  j  ptq-K  J  —  HI  — 

where  v  is  the  kinematic  viscosity.  In  the  limit  D  ■*  •  (necessary  for 
strict  isotropy), 

E(k)  *  wkU(k),  U(k)  *  (D/^*)2<|u(k) |2>.  (2.2) 


where  (  ^  denotes  ensemble  average  and  the  mean  kinetic  energy  per  unit 
mass  is  E(k)dk.  U(k)  measures  the  intensity  of  excitation  per  mode. 
The  energy  balance  equation  is 


(3/3t  ♦  2vx  )E(k)  =  T(k), 


T(k)  =  i  j  J  T(k,p,q)dpdq,  T(k,p,q)  *  T(k,q,p),  (2.3> 


0  0 


where 


T(k,p,q)  s  2ik  Ira{(2*/|sin(p,q)  |  ) ( D/2 w ) **( k^d j ♦k 

*(u*(k)Uj(p)un(^)>)  (k=p+q,  k  =  |k | ,  p  =  |p|,  q  =  |q|), 

T(k,p,q)  =  0  (if  k,  p,  q  cannot  form  the  sides  of  a  triangle).  (2.4) 


To  obtain  (2.4)  from  (2.1),  use  -*  (D/2*)2; Jd‘p  (D  ■*  •)  and  JJd2p  * 
2«//lsin(p,q) |”1dpdq. 

Detailed  conservation  of  energy  and  squared  vorticity  for  each  triad 
interaction  is  expressed  by 

T(k,p,q)  ♦  * ( p , q y *  ^  ^  *  0 y 


k2T(k,p,q)  ♦  p2T(p,q,k)  ♦  q2T(q,k,p)  =  0, 


(2.5) 
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which  can  be  verified  from  (2.4)  with  the  use  of  incompressibility  and 
plane-tria  .gle  identities.  The  overall  conservation  laws 

m  m 

I  T(k)dk  *0,  /  k2T(k)dk  *  0  (2.6) 

0  0 

follow  from  (2.5).  Conversely,  (2.6)  implies  (2.5),  since  (2.6)  must  hold 
fer  states  in  which  only  a  single  triad  of  wavevectors  have  nonzero  ampli¬ 
tudes  at  an  instant  so  that  the  instantaneous  T(  )  vanishes  outside  the 
triad.  By  (2.5), 

T(p,q,k)/T(k,p,q)  =  (q2-k2)/(p2-q2 ) ,  T(q,k ,p)/T(p,q,k) 

*  (k2-p2)/(q2-k2),  T(k,p,q)/T(q,k,p)  *  (p2-q2)/(k2-p2)  (2.7) 

so  that  only  one  of  the  T(  ,  ,  )  associated  with  a  given  triad  interaction 
is  linearly  independent. 

The  mean  rate  of  transfer  of  kinetic  energy  per  unit  mass  from  wave- 
numbers  below  k  to  those  above  is  I!(k)  =  T(k')dk'.  By  (2.5), 

•  k  k  k  -  • 

n(k)  =  i  J  dk'J  /  T(k',p,q)dpdq  -  i  /  dk'J  J  T(k • ,p,q)dpdq.  (2.8) 
k  0  0  ^  0  k  k 

Th«  first  term  on  the  right-hand  side  is  the  total  rate  of  gain  in  the 

range  k'  >  k  due  to  triad  interactions  with  p,q  <  k,  while  the  second 

term  is  the  total  rate  of  loss  in  the  *ange  k'  <  k  due  to  triads  with 

p,q  >  k.  These  two  classes  of  triad  interactions  are  mutually  exclusive 

and  exhaust  the  interactions  which  contribute  to  net  energy  transfer 

across  k.  Similarly,  the  mean  rate  of  transfer  of  squared  vorticity 

from  below  k  to  above  k  is 

"  k  k  dpdq  k  •  • 

Z(k)  *  /  (k'Tdk'J  /  T(k',p,q)  -  /  (k*)  <ik'J  /  T(k*  ,p,q)dpdq. 
k  0  0  *  0  k  k 


(2.9) 
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Assume  that  the  double  and  triple  movents  at  the  instant  considered 
satisfy  the  similarity  laws 

E(ak)/E(k)  =  a'n,  T(ak,ap,aq)/T(k,p,q)  *  (2.10) 


where  a  is  an  arbitrary  scaling  factor  and  n  is  so  far  undetermined.  The 
scaling  of  T(k,p,q)  in  (2.10)  is  the  same  as  that  of  [E(k)]3^2k~172 
(which  has  the  same  dimensions)  and  corresponds  to  a-independence  of  the 
appropriately  defined  triple-correlation  coefficients  of  the  distribution 

of  the  Fourier  amplitudes  in  the  neighborhoods  of  the  wavenumber  arguments. 

f  k  k  k  o 

Note  that JQ  dp/Q  dq  is  equivalent  to  2 fQ  dp/^dq  in  the  first  term  on 

the  right-hand  side  of  (2.8)  because  or  the  symmetry  of  T(k’,p,q),  while 

Jk*dpJk  dq  in  the  second  tern  is  equivalent  to  2/^  dp/^  dq.  Set  p  =  k/u, 

q  =  ov,  k'  *  pw  in  the  first  terra  and  p  *  k/u,  k'  =  pv,  q  =  pw  in  the 

second  term.  Note  that  du/^  dw  is  equivalent  to  dw/^du  in  the 

first  term  and  J^du/^dv  is  equivalent  to  J^dvJ^du  in  the  second  term. 

Use  (2.10)  with  a  a  k/u,  and  finally  use  (2.7)  to  obtain 

n(k)  a  k(5’3n>/?J  dv/  dw  W  (v,w,n)T(l,v,w),  (2.11) 

0  1  1 


where 


W,(v,w,n)  *  -(w2-v2)*1[(l-v2)/  u(3n‘7)/2du-(w2-l)J  u(  3n_7)/2du].  (2.12) 


Repeating  the  procedure  for  2(k)  gives 


1  - 


Z(k)  =  2k(9-3n)/2J  dv/  dw  W  ( v,w,n)T(l,v,w), 


(2.13) 


0  1 


where 

„  ,  x  ,  2  2.-lr/.  2.  2rW  (3n-ll)/2.  ,  2  2rX  (3n-ll)/2. 

vMv»w*n)  s  -(«  ~v  )  C(l-v  )w  J  u  du-( w  -l)v  J  u  duj. 

2  1  v 


(2.14) 
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Equations  (2.11)  and  (2.13)  er  f su  Il(k)  and  Z(k)  as  integrals  over 
contributions  from  all  the  possible  shapes  of  the  triangles  formed  by 
k',  p,  q  in  (2.8)  and  (2.9).  Since  v  <  1  and  w  >  1,  each  pair  of  values 
vf  w  corresponds  uniquely  to  a  particular  triangle  shape,  by  definition, 
T(l,v,w)  is  zero  if  1,  v,  w  cannot  for®  a  triangle.  The  W  factors  give 
the  weights  of  the  contributions  of  the  different  triangle  shapes  and 
arise  fro®  integration  over  triangle  size. 

If  n  :  5/3,  (2.11)  says  that  n(k)  has  a  value  c  which  is  independent 
of  k.  If  n  =  3,  (2.13)  gives  7.(k)  a  value  n  which  is  independent  of  k. 

By  (2.12)  and  (2.14), 

W/v.w.5/3)  *  -(w2-v?)‘1[(l-v2)ln(w)+(w2-l)ln(v)j, 

«2(v,w,5/3)  =  0,  ^j(v,w,3)  -  o, 

W?(v,w,3)  =  -(w?-v‘)‘1[(l-v2)w2ln(w)+(w2-l)v2ln(v)].  (2.15) 

Thus,  for  each  triangle  shape  individually,  an  n  =  5/3  similarity  range 
yields  a  k-independent  energy  cascade  and  identically-zero  vorticity 
cascade,  while  an  n  *  3  similarity  range  yields  a  k-independent  squared- 
vorticity  cascade  and  identically-zero  energy  cascade. 

The  scaling  of  E(k)  ha3  not  been  used  in  obtaining  these  results. 
Therefore  they  hold  al30  for  mote  general  similarity  ranges  in  which  n 
is  replaced  by  n'  /  .i  in  the  first  equation  of  (2.10). 
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3.  CASCADE  DIRECTIONS 

T(i,v,w)  in  (2.11)  and  (2.13)  represents  a  (signed)  flow  into  the 
middle  wavenumber  cf  the  triad  l,v,w  since  v  <  1  and  w  >  1.  It  is  shown 
in  the  Appendix  that  >*^(^*,5/3)  >  0  and  W^(v,w,3)  <  0.  This  means  that 
the  contribution  of  each  triangle  shape  to  e  in  the  -5/3  range  has  the 
same  sign  as  the  flow  of  excitation  into  the  middle  wavenumber,  while  the 
contribution  to  n  in  the  -3  range  has  the  opposite  sign.  There  is  nothing 
in  the  conservation  properties  by  themselves  to  determine  the  sign  of 
T(l,v,w).  Indeed,  if  a  similarity  range  with  a  given  sign  of  T(l,v,w) 
exists  at  an  instant,  then  a  range  with  the  opposite  sign  is  produced  by 
reversing  the  velocity  everywhere  in  space. 

Physical  interest  attaches  not  to  hypothetical  instantaneous  similar¬ 
ity  ranges  but  to  the  possibility  of  quasisteady  ranges  which  develop 
under  the  dynamical  equations.  In  three  dimensions,  a  hint  as  to  the 
direction  of  cascade  in  the  -5/3  range  comes  from  considering  the  absolute 
statistical  equilibrium  which  would  obtain  if  viscosity  were  zero  and  the 
system  were  truncated  by  removing  all  degrees  of  freedco  with  k  greater 


than  seme  cutoff  wavenumber  k  from  the  dynamical  equations.  The  total 

max 

kinetic  energy  per  unit  mass  is  ^  J.  |u(k)|2,  and  consequently  the  energy 

2  £  -  ~ 

spectrum  in  the  hypothetical  absolute  equilibrium  would  have  the  equiparti- 

2 

tion  fora  U(k)  *  constant,  or  E(k)  ■  k  .  The  -5/3  spectrum  means  that 
higher  wavenumbers  in  the  inertial  range  are  far  below  absolute  equilibrium 
with  lower  wavenumbers  and  it  i9  plausible  that  the  dynamical  interaction 
should  act  toward  producing  equilibrium,  a  state  which  never  can  be  reached 
because  the  viscous  dissipation  provides  a  high-k  sink. 
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In  two  dimensions,  the  absolute  equilibrium  has  a  more  complicated 
structure  because  there  are  two  linearly  independent  quadratic  constants 
of  motion.  The  general  fcrm  of  the  equilibrium  spectrum  is 

U(k)  *  l/(Bk2  ♦  a),  (3.1) 

where  6  and  a  are  constants.  This  is  an  equipartition  distribution*0  for 

2  2 

the  constant  of  motion  £k(0k  ♦  a)|u(k)|  .  The  corresponding  vorticity 

3  ~ 

spectrum  2*k  U(k)  increases  aonotonically  with  k  so  that  most  of  the 
vorticity  in  equilibrium  is  at  wavenumbers  k  .  Since  k  can  be 
arbitrarily  high,  this  suggests  that  a  tendency  toward  equilibrium  in 
an  actual  physical  flow  should  involve  an  upward  flow  of  vorticity  and 
therefore,  by  the  conservation  laws,  a  downward  flow  of  energy.  Thus  if 
the  nonlinear  interaction  does  act  toward  producing  equilibrium,  T(l,v,w) 
should  be  typically  negative.  A  sinpler  and  cruder  statement  is  that 
T(l,v,w)  should  be  negative  because  that  represents  a  statistically 
plaur'ble  spreading  of  the  excitation  in  wavenumber:  out  of  the  middle 
wavenumber  into  the  extremes. 

Supporting  evidence  is  provided  by  the  initial  growth  of  energy 
transfer  in  turbulence  whose  initial  statistical  distribution  is  Gaussian. 

The  exact  expression  for  this  in  two  dimensions  has  been  obtained  by 

8  9 

Reid  and  Ogura.  It  is 

[dT(k,p,q)/dt]0  r2»2k2dkpqC2akpqU(p)U(q)-bkpqU(q)U(k)-bkqptJ(k)U(p)]. 


(3.2) 
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Here 

b.  s  2pk*i(xy-z>*-2z3),  2a.  =  b.  «-b.  t 

fcpq  3  *  »<pq  *pq  *qp 

2  1/2 

d,  =  k/(l-x  )  =  diameter  of  circumscribed  circle.  (3.3) 

kpq 

wnere  x,  y,  z  are  the  interior  angles  opposite  the  triangle  sides  k,  p,  q. 
The  coefficients  obey  the  identities 

a,  >  0,  d,  *  d  .*  d  .  ,  k2b.  *  p2b  l.  » 

%pq  kpq  pqk  qkp*  kpq  F  pkq* 

2k2a.  *  p2b.  ♦  q2b.  (3.4) 

kpq  *  kpq  *  kpq 

whence 

=  «12->‘J>/(<12-P2>.  V'V)  *  2(p2-k2)/(p2-<l2).  0.5) 

Now  suppose 

U(p)  =  (p/k)"rU(k),  U(q)  =  (q/k)“rU(k).  (3.6) 

Equations  (3.2)-(3.5)  yield 

CdT(kiPfq)/dt]0  =  2*2k2dkp4akpq(pq/k2)"r[U(k)]2 

*  (l-vr(w2-l)/(w2-v2)  -  wr(l-v2)/(w2-v2)},  (3.7) 

where  v  =  pA  ,  w  s  qA.  It  is  shewn  in  the  Appendix  that  when  v  <  1  and 
w  >  1  the  curly  bracket  in  (3.7)  is  >  0  if  0  <  r  <  2  and  <  0  if  r  <  0  or 
r  >  2.  Thus  the  initial  growth  of  T(k,p,q)  gives  a  positive  flow  into  the 
middle  wavenumber  k  if  r  falls  between  the  limits  (0,2)  set  by  the  extreme 
absolute  equilibrium  distributions  o  =  0  and  a  =  0  in  (3.1).  If  r  is  out- 
sice  these  limits,  there  is  net  flow  out  of  k.  The  -5/3  and  -3  similarity 
ranges  correspond  to  r  =  8/3  and  r  =  4,  both  of  which  yield  net  outflow. 
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4.  CONJECTURES  ON  QUASI STEADY  STATES 

Are  the  formal  -5/3  and  -3  similarity  ranges  asymptotic  limits  of 

states  which  can  arise  physically?  Suppose  that  an  infinite  fluid  is 

excited  by  isotropic  stirring  forces  confined  to  k  k.,  where  k^  is  a 

characteristic  input  wavenumber.  Let  the  stirring  forces  supply  energy 

2 

at  a  steady  rate  c  and  squared  vorticity  at  a  steady  rate  n  'v*  2k?  c. 

More  general  ratios  n/c  are  interesting  but  will  not  be  considered  here. 

The  preceding  analysis  suggests  that  if  the  input  continues  for  a  suffic- 

iently  long  tine  and  the  Reynolds  number  [E(k^)A^]  /v  is  large  enough, 

a  quasisteady  state  may  be  set  up  in  which  an  approximate  -3  vorticity- 

transfer  range  carries  most  of  the  squared-vorticity  input  up  to  k  >>  k^, 

where  it  is  dissipated  by  viscosity,  while  an  approximate  -5/3  energy- 

transfer  range  carries  most  of  the  energy  input  down  toward  zero  wave- 

number.  The  -5/3  range  can  be  only  quasisteady  because  its  lower  end 

keeps  moving  down  to  ever-lower  wavenumbers,  a  wavenumber  k  <<  k^  being 

2  -1/3 

reached  in  a  tine  t  ^  (ck  )  according  to  energy  conservation.  As 
t  •*  «*,  the  rate  of  transfer  of  squared  vorticity  from  k  t  k.  to  lower 
wavenumbers  decreases  steadily  toward  zero  and  the  energy- transfer  range 
approaches  the  asymptotic  -5/3  dynamics  ever  more  closely. 

The  formal  similarity  ranges  can  represent  asymptotic  quasisteady 
states  only  if  the  cascades  are  sufficiently  local  in  wavenumber.  The 
questions  involved  here  are  the  same  in  two  and  three  dimensions.  Local 
transfer  in  the  -5/3  range  is  plausible  for  the  reasons  given  by  Kolmo¬ 
gorov.  The  transfer  is  associated  with  the  distortion  of  the  velocity 
field  by  its  own  shear.  The  cascade  !!(k)  through  a  given  wavenumber  k 
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in  a  -5/3  range  is  expected  to  be  negligibly  affected  by  wavenumbers  <<  k 
because  the  integral  fQ  k'L(k)dk,  which  measures  the  mean-square  shear, 
converges  at  k  =  0.  It  is  expected  to  be  negligibly  affected  by  wavenumbers 
>>  k  because  the  vorticity  associated  with  those  wavenumbers  fluctuates 
rapidly  in  space  and  time  and  gives  an  effective  shear  across  distances 
of  order  k~*  which  is  small  ccmpared  to  the  shear  associated  with  the 
wavenumbers  ^  k. 

On  the  other  hand,  the  s 4uared-vorticity  spectrum  in  the  -3  range  is 
•  k-i  so  that  each  octave  below  a  given  wavenumber  k  contributes  the  sane 
amount  to  the  mean-square  shear  and  the  latter  diverges  logarithmically 
toward  small  k.  This  means  that  transfer  in  the  -3  range  is  not  local 
when  a'l  triad  Interactions  are  admitted.  It  would  seem  plausible  that 
when  this  range  occurs  as  a  quasisteady  state  the  power  law  is  modified 

m 

by  logarith^c  corrections.  Subject  to  such  corrections,  the  -3  range  can 

3  1/6 

be  expected  to  extend  up  to  k  “v  k^  s  (n/v  )  ,  at  v.iich  wavenumber  the 

integrated  rate  of  dissipation  of  squared  vorticity  by  viscosity  reaches 
the  order  of  n.  At  higher  wavenumbers,  E(k)  is  expected  to  fall  off  at 
a  faster-than-algebraic  rate. 

2 

The  corresponding  viscous  dissipation  of  kinetic  energy  is  'x.  nAj  , 

2 

which  implies  c(k.  Ad>  .  Thus  c^/c  ■*  0  if  v  ■»  0,  in  marked  contrast 

to  three-dimensional  turbulence  where  the  energy  loss  becomes  independent 
of  v  as  v  ■*  0,  These  considerations  imply  that  the  entire  energy  input 
c  is  carried  down  toward  k  =  0  without  viscous  loss  in  the  limit  of  infinite 
Reynolds  number.  The  result  is  not  directly  applicable  to  meteorological 
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flows  because  the  constraints  which  render  the  latter  two-dimensional 
break  down  at  sufficiently  high  k. 

The  actual  correctness  of  the  strict  -5/3  asymptotic  inertial- 
range  law  is  not  established  beyond  doubt  in  three-dimensional  flow, 
and  the  arguments  for  similarity  ranges  in  two  dimensions  are  substantially 
less  secure.  The  present  paper  has  demonstrated  some  elementary  consistency 
properties,  but  this  does  not  show  that  the  similarity  ranges  actually 
exist. 

One  important  difference  between  two  and  three  dimensions  is  the 

existence  of  an  infinite  number  of  local  inviscid  constants  of  motion 

in  the  former:  the  vorticity  of  each  fluid  element.  This  implies  that 

inertial  forces  alone  cannot  produce  universal  statistical  distributions 

in  the  similarity  ranges,  independent  or  the  statistical  distribution  of 

the  driving  forces.  In  three  dimensions  there  are  also  an  infinite 

number  of  inviscid  constants  of  motion:  the  circulations  about  all  closed 

curves  moving  with  the  fluid.  However  a  given  closed  curve  is  expected 

to  stretch  and  migrate  in  complicated  fashion  through  the  fluid  with  the 

passage  of  time  so  that  it  is  reasonable  to  expect  that  the  circulation 

invariance  does  not  impose  effective  constraints  on  n-variate  distribution 

(the  latter  corrected  by  a  logarithmic-type  function  of  k/k.) 
functions  for  small  n.  If  (1.1)  and  (1.?)  are  realized  in  two  dimensions^, 

it  is  tc  be  expected  that  C  and  C*  are  not  universal  constants  but  depend 

on  the  character  of  the  driving  force?,  whatever  may  be  the  situation  in 

three  dimensions.  A  further  point  is  that  the  nonlocalness  of  transfer  in 

in  itself 

the  -3  range  suggestsAthat  cascade  there  is  not  accompanied  by  degradation 
of  higher  statistics  in  the  fashion  usually  assumed  in  a  three-dimensional 
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Kolraogorov  cascade.  This  is  consistent  with  a  picture  of  the  transfer 
process  as  a  clumping-together  and  coalescence  of  siirilarly  signed  vortices 
with  the  high  wavenumber  excitation  confined  principally  to  thin  and  in¬ 
frequent  shear  layers  attached  to  the  ever- larger  eddies  thus  formed. 

In  connection  with  the  sign  of  c  in  the  -5/3  range,  it  should  be 

noted  that  a  positive-c  range  extending  from  the  input  wavenumbers  up 

to  the  dissipation  range  would  be  physically  unrealizable.  The  viscous 

dissipation  would  remove  squared  vorticity  and  kinetic  energy  in  a  ratio 

nd/Cd  which  would  greatly  exceod  the  ratio  at  which  these  quantities 

were  cascaded,  since  the  letter  ratio  goes  to  zero  as  the  asymptotic 

-5/3  structure  is  approached.  Thus  an  upward-transfering  -5/3  range 

could  not  exist  as  a  quasisteady  state.  The  inconsistency  of  such  a 

2 

range  has  previously  been  demonstrated  by  Lee. 

Suppose  now  that  the  fluid  is  confined  to  a  finite  domain  and  that 
♦he  lowest  wavenumber  allowed  by  the  boundary  conditions  is  kQ  <<  k.. 

The  conjecture  is  offered  here  that  after  the  -5/3  range  reaches  down 
to  wavenumbers  *  kQ  the  downward  cascade  froei  k^  continues  and  the  energy 
delivered  to  the  bottom  of  the  range  pil ?s  up  in  the  mode  kQ.  As  the 
energy  in  kQ  rises  sufficiently,  modification  of  the  -5/3  range  toward 
absolute  equilibrium  is  expected,  starting  at  the  bottom  and  working  up 
to  progressively  larger  wavenumbers. 

Come  support  for  the  idea  of  energy  piling  up  in  kQ  comes  from  con¬ 
sidering  the  absolute  equilibrium  ensembles  (3.1).  Suppose  that  the 
wavenumber  range  is  truncated  from  below  at  kQ  ana  from  above  at  a  wave¬ 
number  k  .  Let  the  mean  energy  and  mean-square  vorticity  have  specified 
ndx 
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2 

values  E  and  Q  =  2k  ^  E.  The  ratio  and  signs  of  8  and  a  in  (3.1)  depend 

on  the  relative  values  of  krt,  k,,  and  k  .  In  particular, 

0’  1*  max  r  ' 


a  *  0  if  k, 


(4.1) 


6s0if  ki2s7(kL*ko2)- 


(4.2) 


2  2 
For  values  of  k^  between  (4.1)  and  (4.2),  a  >  0,  8  >  0.  For  k^  less  than 

(4.1),  8  >  0,  -8kQ2  <  a  <  0.  If  kL-k  <<  kQ,  then  8kQ2  ♦  a  «  8k2  and  E(k) 

2  2 

has  a  sharp  peak  at  k  =  kQ.  For  k^  greater  than  (4.2),  a>0,  -a/k(nax  <  8 

2 

<  0,  and  k  E(k)  shows  a  sharp  peak  at  k  =  k  if  k  -k,  <<  k  .  The 

r  r  max  max  1  max 

values  k  <  k  and  k.  >  k  are  impossible.  These  results  are  all  for 
xu  x  rr*x 

a  continuous  spectrum  of  allowed  wavenumbers.  When  the  discreteness 
associated  with  a  finite  fluid  is  taken  into  account,  the  sharp  peak  in 
E(k)  is  modified  so  that  the  lowest  mode  kQ  singlehandedly  carries  most 
of  the  total  kinetic  energy  if  k^-kQ  is  much  less  than  the  mode  separation  8k 
of  the  lowlying  degrees  of  freedom. 

These  results  suggest  that  a  piling  up  of  energy  in  kQ  under  a 
steady  input  would  represent  a  plausible  way  for  the  wavenumbers  <  k^  to 
seek  an  absolute  statistical  equilibrium  of  the  kind  that  corresponds  to 
very  large  E/Q.  The  phenomenon  is  analogous  to  the  Einstein-Bose  conden¬ 
sation  of  a  two-dimensional  quantum  gas. 
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APPENDIX.  PROOFS  OF  INEQUALITIES 

Take  w>l,0<v«l  throughout  this  Appendix.  Let 
I(v,w)  *  (l-v2)ln(w)+(w2-l)ln(v). 

Then  I(l,w)  =  0,  3I(v,w)/3v  =  (w2-l)/v  -  2vln(w),  and 
[dl(v,w)/>v]v=1  s  (w2-l)  -2  ln( w)  i  F(w) . 

Now  F(l)  =  0,  dFlw)/dw  =  2w  -  2/w  >  0.  Therefore  w2-l  >  2  ln(w)  and, 
since  v"1  >  v,  it  follows  that  3l(v,w)/3v  >  0,  so  that  I(v,v)  <  0.  This 
establishes  W^(v,w,5/3)  >  0. 

Let 

J(v,w)  *  (l-v2)w2ln(w)*(w2-l)v2ln(v). 

Then  J(l,w)  =  0,  3J(v,w)/3v  *  (w2-l)v[lt21n(v) ]-2vw^ln(w) .  Now2w2ln(w) 

2 

>  w  -1  and  v  -  vt2vln(v),  the  first  inequality  readily  following  uoon 
differentiation  and  the  second  foil*"  .xng  from  ln(v)  <  0.  Therefore, 

3J( v,w)/3v  <  0,  J(v,w)  >  0,  whence  V?2(v,w,3)  <  0. 

Let 

K(v,w,r)  =  l-vr(w2-l)/(w2-v2)  -  wr(l-v2)/(w2-v2). 

r?  -1  r?  25 

3K(v,w,r)/3r  =  [v  (w  -l)ln(v  j  -  w  ll-v  )ln(w) ]/(*.•  -v^) , 

considered  as  a  function  of  r,  is  the  sum  of  a  positive  term  of 

raondtonically  decreasing  magnitude  and  a  negati .e  term  of  nonotonically 

increasing  magnitude.  Therefore  it  has  just  one  zero,  and  it  is  positive 

is  r  ■»  negative  as  r  ♦  Since  K(v,w,r)  has  zeroes  at  r  =  0,  r  =  2, 

it  follows  that  K(v,w,r)  is  >  0  for  0  <  r  <  2  and  <  0  for  r  <  0  and  r  >  2. 

This  establishes  the  sign  of  the  right-hand  side  of  (3.7). 
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FIGURE  CAPTION 

FIG.  1.  Part  (a)  represents  the  cascade  of  oxcitation  through  the 
wavenumber  spectrum  by  means  of  elementary  triad  interactions.  Part  (b) 
represents  a  pair-interaction  cascade,  which  is  a  valid  simplification  in 
three  dimensions  but  not  in  two  dimensions. 
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II  ABbTMAfT 


Two-dimensicnal  turbulence  has  both  kinetic  energy  and  mean-square  vorticity 
as  inviscid  constants  of  motion.  Consequently  it  admits  two  formal  inertial 

ran*'s  Mk)  a.  cJ/V5/3  and  £(k)  a.  nJ/V\ 

where  c  is  the  rate  of  cascade  of  kinetic  energy  per  unit  mass,  n  is  the  rate 
of  cascade  of  mean-square  vorticity,  and  the  kinetic  energy  per  unit  mass  is 

J0*E(k)dk. 

The  -S/3  range  is  found  to  entail  backward  energy  cascade,  from  higher  to 
lower  wavenumbers  k,  together  with  zero  vorticity  flow.  The  -3  range  gives 
an  upward  vorticity  flow  and  zero  energy  flow.  The  paradox  in  these  results 
is  resolved  by  the  irreduoioly  triangular  nature  of  the  elementary  wavenumber- 
interactions.  ihe  formal  -3  range  gives  a  nonlocal  cascade  and  consequently 
must  be  modified  by  logarithmic  factors.  If  energy  is  fed  in  at  a  constant 
rate  to  a  band  of  wavenumbers  ^  k^  and  the  Reynolds  number  is  large,  it  is 
conjectured  that  a  quasisteady  state  results  with  a  -S/3  range  for  k  <<  k^ 
and  a  -3  range  for  k  >>  k^,  up  to  the  viscous  cut-off.  The  total  kinetic 
*\.ergy  increases  steadily  as  the  -5/3  range  pushes  to  ever-lower  k,  until 
scales  the  size  of  the  entire  fluid  are  strongly  excited.  The  rate  of 
energy  dissipation  by  viscosity  decreases  to  zero  if  kinematic  viscosity  is 
decreased  to  zero  with  other  parameters  unchanged. 
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